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I describe a method for determining the coefficients of scalar integrals for one-loop amplitudes in quantum 
field theory. The method is based upon generalized unitarity and the behavior of amplitudes when the free 
parameters of the cut momenta approach infinity. The method works for arbitrary masses of both external and 
internal legs of the amplitudes. It therefore applies not only to QCD but also to the Electroweak theory and to 
quantum field theory in general. 



I. INTRODUCTION 

For many years, high energy physics has been looking for- 
ward to the beginning of the LHC program. There are high 
expectations that experiment will answer a number of out- 
standing theoretical issues, such as the nature of electroweak 
symmetry breaking and the resolution of the hierarchy prob- 
lem. However, the success of the experimental program will 
depend upon the development of better theoretical tools. Our 
ability to identify new physics at the TeV scale will depend, 
in large part, on our ability to accurately describe and separate 
out the known physics of the Standard Model. 

It is well established that at the Tevatron and other col- 
liders that leading order calculations of QCD processes are 
insufficient for an accurate description of hadronic interac- 
tions. Next-to-leading order (NLO) corrections in QCD are 
often quite large. The corrections typically change both the 
predicted magnitude of the scattering process and the shapes 
of distributions. 

At the Tevatron, which has little reach into the multi- 
hundred GeV range, it has been found that computing NLO 
QCD corrections is generally sufficient to reduce the theoret- 
ical uncertainty of a calculation to about 10%. At the LHC, 
however, with its reach extending above the TeV scale, it has 
been found that electroweak corrections can become quite 
large. It is therefore becoming important to develop tools 
for performing higher-order corrections in the full SU(3) (g> 
ST/ (2) (gi U(l) Standard Model. 

In the past, many of the most exotic computational methods 
focused on QCD not just because QCD corrections are large 
even at relatively low energies, but also because working with 
massless quarks and gluons simplifies expressions sufficiently 
to allow complicated algorithms to be worked out. Recent ad- 
vances now provide a framework in which massive theories 
can access the sophisticated advances developed for QCD cal- 
culations. 

The workhorse in perturbative calculations has long been 
the Feynman diagram approach, which approach systemat- 
ically includes all perturbative effects. However, the rapid 
growth in the number of Feynman diagrams with the num- 
ber of external legs, even for tree-level amplitudes, limits the 
practical application of this technique. In addition, Feynman 
diagram calculations are subject to large cancellations among 



individual diagrams, especially when working with gauge the- 
ories. 

One tool that has long been used in QCD calculations to ob- 
tain compact results and largely avoid unphysical singularities 
and gauge cancellations is the helicity method HI 12, El IH 0, 
la 13, |Sh • When combined with recursion relations, both off- 
shell IIS and, more recently, on-sheU HI QJ, GJ, H the 
helicity method can be used to generate compact expressions 
for for very complicated scattering processes. 

Another very powerful technique that has greatly simpli- 
fied the calculation of loop amplitudes in QCD is the unitarity 
method EHILlllltllElEIlEi. An essential fea- 
ture of the unitarity method is that it sews together tree-level 
amplitudes into loop amplitudes. Thus, efficient techniques 
for computing tree-level amplitudes which give compact ex- 
pressions, like the use of recursion relations and the helicity 
method, directly benefit the unitarity method of computing 
loop amplitudes. 

There have been a number of new developments that have 
significantly enhanced the power of recursive methods and 
the unitarity method. These include the use of maximally 
helicity- violating (M HV) vertices in recursion relations at 
tree-level SIMS EE, E3 El El El El El and in 
loops | 34l35Ll36Ll37ll . Unitarity methods have been improved 
by the use of the holomorphic anomaly [38, 39), 4C| to evalu- 
ate cuts and the use of complex momenta M4ll ,l42l ,l43l 144 14511 
within the context of generalized unitarity II46L 1471 148L I49L l5(X 
|5jJ,|52j], which allows for the use of multiple cuts in a single 
amplitude. 

These techniques can be combined into a "unitarity boot- 
strap" 15 3l |54J], a systematic, recursive approach to making 
high-multiplicity QCD calculations practical. The bootstrap 
combines the use of four-dimensional unitarity to determine 
the logarithmic and polylogarithmic terms in the loop ampli- 
tude with on-shell recursion relations to determine the rational 
contributions. 

Recently, Forde lf55ll has described an efficient method for 
extracting the coefficients of the triangle and bubble func- 
tions in amplitudes where only massless particles circulate 
in the loops. Combined with the result of Ref. l4ltl . this is 
sufficient to determine the cut-constructable part of one-loop 
QCD amplitudes. In this paper, I will extend Forde's for- 
malism and define methods to extract the coefficients of all 
loop-integral functions in massive as well as massless theo- 
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ries. This will then allow this method to be used in the full 
SU(3)®SU(2)®U(1) Standard Model and even in the Su- 
persymmetric Standard Model. 

The plan of the paper is as follows: In Section (HJl, I will 
discuss the notation used and the application of spinor helic- 
ity conventions to massive complex momenta. In Section dlTll i. 
I will present an overview of the formalism which allows 
the simple extraction of the loop integral functions. In Sec- 
tion (ITVD . I will discuss the extraction of the box functions; 
in Section (0 the extraction of the triangle functions; in Sec- 
tion ( 1VII ), the bubble functions; and in Section (TVIQ . the tad- 
poles. Finally, I will comment on my results and draw my 
conclusions. 



n. NOTATION AND CONVENTIONS 

An important feature of the method described in this paper 
is that it involves on-shell amplitudes where some of the on- 
shell momenta are complex. Therefore, I will describe the 
conventions used to handle complex momenta and, especially, 
the spinorial representations of complex momenta. 



A. Spinor Representations of Real Momenta 

I treat real momenta using the spinor-helicity formalism of 
Mangano and Parke [8], which I will summarize briefly. Let 
be a massless momentum in Minkowski space, and y(p) 
the Dirac spinor representing a massless fermion of momen- 
tum p, 



Some important identities are: 



PY(p) = 0, P 2 = 0. 



The two helicity states of y(p) are given by 



Y±(P)= ^¥(p)(l±7s) = ■ 



The phase convention is chosen such that 



\p ± ) = \p t ) c , (p ± \= c (p^ 



(i) 



(2) 



(3) 



where the c indicates charge conjugation. To render expres- 
sions more compact, I adopt the notation 



{pq) = (p \q + ) , 
[pq] = (p + \q~) ■ 



(4) 



(p + \q + ) 
(pp) 
(pq) 
(pq) [qp] 
(p ± \f\p ± ) 
P 

7,(p ± \f\q ± ) 
(AB) (CD) 
[AB][CD] 



(p-\q-)= 
\PP]= o, 

-(qp) 

2p-q, 

\P + )(P 



[pq] = - [qp] 



(5) 



= 2 { P' 



z )(q 



(AC) (BD) 
[AC][BD]- 



^ \P )(P I > 

1+ I^X/'I 

- (AD) (CB) , 
[AD] [CB] . 



The last two identities are known as the Schouten identities. 

This formalism can alternatively be expressed in terms of 
Weyl-van der Waarden (WvdW) spinors [56]. Indeed there is 
a one-to-one correspondence between the spinors defined here 
and WvdW spinors: 



(<Px) 



MP 



<t>A 




X A 

■dW i 



[%<!>} 

o 

a H,AB 



MP 



(0 A o), 
(o x A ), 

* ($X)wvdW ' 



(6) 



AB 





Where the MP and WvoW subscripts denote the "Mangano & 
Parke" and "Weyl-van der Waerden" notation for spinor prod- 
ucts, respectively. 



1. Fermion Wave Functions for Real Massless Momenta 

These are essentially trivial, since the notation is defined in 
terms of massless Dirac spinors. The only point of convention 
is that, for the sake of defining helicity, all external particles 
are taken to be outgoing. This means that we identify 

|p ± )=v±( /3 ), (p ± \=u ± (p). (7) 

The phase convention defined in Eq. (O means that 

u ± (p) =v T {p), u±{p)=v T (p). (8) 



2. Massless Vector Boson Wave-functions for Real Momenta 

Massless spin-1 particles have two physical polarization 
states. The standard practice in spinor helicity methods is to 
use the light-like axial gauge, in which the polarization vec- 
tors are defined in terms of both the momentum vector k and 
a reference vector g. The gauge invariance of the spin-1 field 
manifests itself in the arbitrariness of the reference momen- 
tum. For an outgoing vector field of momentum k, the helicity 
states are can be written as 



(k+ \-f 



V2(gk) 



(k- \f 



V2[kg] 



(9) 
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These polarization vectors have the usual properties 



e ± -e ± 



o, 



(10) 



V = — CT^ V 



The arbitrariness of the choice of g can be seen by examining 
the difference between two choices of g. Consider contracting 
a polarization vector with some random vector, 

(e + "(*;*i) -e+»(k;g 2 ))(a-\Yn\b-) 



V2[kb] 
V2[kb] 



V (gi k ) (82k) 

(agi) (gik) - (ag 2 ) (gik) 



If 



(gik) (gik) 
(ak) (g 2 gi) 
(gik) (gik) 
Hgi) 



(ID 



V2( gl k) (g 2 k) 



(a \Yn\b ) 



Thus, the difference between the polarization vectors gener- 
ated by two choices of g is proportional to W 1 . Since gauge 
bosons couple to conserved currents, this is a pure gauge term 
and does not contribute to the amplitude. 



B. Massive Real Momenta 



Note that the different helicity states are not related by com- 
plex conjugation, as they are in the case of real momenta. As 
in the case of real momenta, however, these spinor representa- 
tions can be used directly as massless fermion wave function 
and in massless vector wave-functions (see Section dll Al l). 



D. Massive Complex Momenta 

Let us assume that we have an on-shell complex momentum 
1 with mass m parametrized in terms of real massless four- 
momenta an d 



wy 



2t At% ■ y 



(v-\y"\x-)- 



(14) 



I can trivially decompose £^ into two massless complex mo- 
menta, 



e1=yx» + ™r + ^(x-\f\v-) + ^(v<-\f\x-) , 



l i 



4tx-V 

and then find spinor representations of £^ and > 



(15) 



A massive real momentum can be represented as the sum 
of two massless momenta. There is great freedom in choos- 
ing this decomposition, which leads to a variety of choices for 
representing massive fermion and vector wave functions. The 
construction of helicity basis wave functions for massive par- 
ticles is discussed thoroughly in Ref. lf56ll and is not repeated 
here. One can even choose the decomposition so that the mas- 
sive wave functions are eigenstates of the helicity projection 
operator. 



C. Massless Complex Momenta 

The massless complex momenta that I will encounter will 
be defined in terms of the spinor representations of two real 
massless momenta x^ an d W* 1 , 

£» = yx»+ w r+ t -(x-\rw-) + ^{v-\f\x-)- 

(12) 

One can define a spinor representation of this complex mo- 
mentum as 



(x + \ + (v" 



(i | = t(x \+w(y- 



\r)= 1\Z~) + \V~), V + )= t\ X + )+w\y+) 









(£1 | - t(x \+w(\jf | 


K> 




-hr> 


K)= t\ X + ) + w\Y+) 


(4\ 


m 1 

~~ lxvV x 




{ l i\ = —, — r<V~ 


V2) 


m , 

" [xv\ |z 


■> 




> 






+ \ti)('i\- 



(16) 



With this parametrization, one immediately finds that 

(kh) = - (hh) = m, [£ih] = - [hh] = m. (17) 



1. Fermion Wave Functions 



If the internal massive particle is a fermion, we must con- 
struct wave functions which obey the Dirac equation. Equa- 
tion (TTTb which implies that 



(18) 



(13) 



l\l+) = -m \q) 


t\q) = 


-m \£+) , 


l\t[) = 


t\4) = 


m \t[) , 


(£+\l=-m(q\ 


(%\*= 


-m{£+\ , 


(£i\t= m(£+\ 


($\* = 


m(q\, 
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which means that the Dirac spinors can be written as: 



2. Massive Gauge Boson Wave Functions 



h(4> = 

lnW> = -|4 + > + K) 

(« T (0|= Wl + <4| 
<v T w| = -(4| + <^| 



h(4X 



K) + I4"X 



(19) 



i I ' 



Note that I label the spin states | / I. rather then ±. This 
indicates that the decomposition of £ into l\ and £ 2 defined 
above does not yield wave functions that are eigenstates of 
the helicity projection operator. Since states with complex 
momenta are necessarily internal states, I do not need helicity 
projections, I only need to sum over the spin states. One can 
easily verify that these spinors obey the Dirac Equation, 



(20) 



the standard normalization conditions, 

{ui{l)\ Uj {l)) = 2mS,j (u i (£)\v j (£))=0 i,je{U} 
(v i (£)\v j (£)) = -2m8 ij (v i (£)\u j (£))=0 

(21) 

and combine to form the standard projection operator, 



/l«t(*)> = 


m | 


l\ui{£)) = 


m \ui{t)) , 


(Ht)) = 


— m |v-f(£)) 


t\*M) = 


-m \vi(£)) , 




m (w| {£) | 




m (ui{l)\ , 


{v^)\l = 


-m(v T (*)| 


<nW|/ = 


-m (v|(f) , 



E !«,(*)> <«,(*)!= £ 



4 4 



|4~X4 I |4 X4~| 



4 4 



4 4 



E |v,w><v,wi = E 

<'e{U} M+-} 

- I4 + X^rl- KX41 + K)(4l+ |4")(^I 
= /-m, 

(22) 

where I can use the Schouten identities to show that 

I4~X41 + KX4I- KX41- \ i i)( i i\= m - ( 23 ) 



- (a£i) (£ 2 b) + (a£ 2 ) (£\b) = - {ab) (£ 2 £i) = m (ab) 
- [a£ 2 ] [Mi] + [a£{\ [£ 2 b] = - [ab] \£ x £ 2 ] = m [ab] . 



(24) 



In this way, I have defined orthogonal spin states for mas- 
sive fermions with complex momenta. Although they are de- 
fined in terms of momentum spinors of definite helicity, these 
spin states are not eigenstates of the helicity projector. 



Massive gauge bosons have three physical spin states. For 
outgoing particles, the wave functions, or polarization vectors, 
can be written as 



e T " = 



e lM = 



{^W 2 ) 
V2[£i£i] 

^«-4) 



(WI4 + ) 

V2m 



(25) 



As with the fermionic wave-functions, I do not use the stan- 
dard symbols of {H — 0} since these spin states are not eigen- 
vectors of the helicity projector. They are, however, orthonor- 
mal and display all of the usual properties expected of polar- 
ization vectors: 







1 £ T £* = 0, 

£ l £^=-l £^£^= £"^ • £* = 0, 

£*£ T = e'-e^ e'e* = -1, 

£»£ v 



(26) 



E. Scalar Loop Integrals 

To establish my sign convention, I will define the scalar 
loop integrals as follows, 

F = (f-mj), 

F 1= ({£ + K l f-m\) 1 

F 2 = ({£ + K 2 f-m 2 2 ), 

F 3 = ((£ + K 3 ) 2 -ml), 

Do{mo\K\,mi\K 2 ,m 2 \K' i ,m^) 

H47rj J (2n) D F Q FiF 2 F?, ' 



C {m ;Ki,mi;K 2 ,m 2 ) = -i(4n 



p/2 



1 



(27r) D F FiF 2 ' 

d D £ 1 



(27) 



Note the sign convention on the momenta in the propaga- 
tors. If, because of routing or some other convention, one 
of the momenta in the denominator is negative, say F\ = 
({£ — K\) 2 — OTj), the loop integral functions would be writ- 
ten as Z?o( w Oj — K\ ,mi), Co(wq, — K\ ,m\,K 2 ,m 2 ), etc. 
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III. METHODS 

In theories of with only massless particles propagating in 
loops, like QCD, it has been shown that any one-loop integral 
can be decomposed into a sum of scalar box, triangle and bub- 
ble loop integral functions that can be constructed from cuts 
using unitarity and a set of rational terms 11611 



\—loop 
n.QCD 



2f 



(4k) 



2-e 



(28) 



where Bq represent the bubble integrals, Cq the triangle 

Ik) 

integrals and D the boxes. The fact that only box func- 
tions are needed is because higher point functions can be 
written as sums of the boxes formed by pinching vertices 
together ||57[ l58tl . The rational terms, 3% n cannot obtained 
from cuts, but can derived from on-shell recursion rela- 
tions OCIBH 

In theories with massive particles, this formula must be aug- 
mented to include tadpole functions, Aq 



, l — loop 



>2e 



(29) 



In standard Passarino-Veltman reduction 15911 . one works in 
D = 4 — 2e dimensions and the coefficients of the loop inte- 
gral functions depend on the dimensional regulator e. Ratio- 
nal terms develop when e-dependent pieces of the coefficients 
multiply poles in e from the loop integral function. In the 
methods described below, the cuts will be evaluated in D = 4 
dimensions, thereby missing those rational terms, which are 
assigned to ffl n . ffl n will then be determined by loop-level on- 
shell recursion relations as part of the unitarity bootstrap. 

The basic idea of the unitarity method is that when one cuts 
a one-loop amplitude, the terms on opposite sides of the cuts 
are tree-level amplitudes. These tree level amplitudes can be 
efficiently computed using helicity methods. However, he- 
licity methods become very messy when they are extended 
into D = 4 — 2e dimensions. It is desirable, therefore, to per- 
form cuts in D = 4 dimensions. Doing so, however, misses 
various rational terms, as described above, which must be ac- 
counted for somehow. In pure QCD, a detailed knowledge of 
the soft and collinear factorization properties of amplitudes is 
sufficient to construct the missing rational terms. In massive 
theories, however, it is possible that there are rational terms 
that cannot be constructed in this way lfl7ll . T he development 
of on-shell recursion techniques 11 ll Il2l Il3l 11411 to compute 
the rational terms allows the reliable use of four dimensional 
unitarity cuts in massive theories as well. 

In generalized unitarity, multiple cuts are made on the 
loop amplitude, dividing it into multple tree amplitudes. The 
method I will describe below uses generalized unitarity to 



determine the coefficients of loop integrals in terms of the 
tree-level amplitudes that appear at the vertices of the cut 
diagrams. Different approaches to determining loop inte- 
gral coefficients include algebraic solution 160116111 ). using D- 
dimensional cuts 162116311 .and using standard unitarity cuts via 
the holomorphic anomaly l64ll . 



IV. FOUR PARTICLE CUTS AND THE EXTRACTION OF 
SCALAR BOX COEFFICIENTS 

The use of quadruple cuts, shown in figure Q] to extract the 
box coefficient was first derived in Ref. [41]. I re-derive the 




FIG. 1 : a quadruple cut. 

result here, using a more efficient momentum parametrization 
and extending the procedure to explicitly permit masses on the 
internal legs. The basic procedure is to cut all four internal 
legs, and put the cut legs on shell, 



{2Tz)8{{l + Ki) 2 -m} 



Imposing the cut conditions, 



f-ml = 0, 
(£ + K 2 ) 2 -m[ 



0. 



-*i) 2 -m? = 0, 
Ki -K A f -m\ = 0, 



(30) 



(31) 



singles out a unique box configuration. Other boxes satisfy a 
different set of cuts and lower point diagrams cannot satisfy 
all four cuts. The contribution of this box topology to the one- 
loop amplitude is obtained by sewing together the four tree- 
level diagrams formed by cutting the loop propagators (See 
figure |2j) 

In D = 4 dimensions, the four delta-function constraints 
completely determine the loop momentum V 1 . There are two 
solutions and,in general, £ is complex valued. The coefficient 
d m , corresponding to scalar box D Q , , is given by averaging 
the product of the four tree-level diagrams, evaluated with the 
solutions to the loop momentum, 



d k 



(32) 



6 



Tk, 





FIG. 2: Tree level amplitudes contributing to a particular box topol- 
ogy- 



solutions for I are: 



at 



S 2 (7i2 -Si) + (712 - S 2 ) ml - Ynm 2 2 + S 2 m\ 
H 2 -SiS 2 

S\ (7i2 - S 2 ) + (7i2 - Si ) ml - Yn m\ + Si m\ 
j\ 2 -SiS 2 



2(ai-l)K\-K 4 + 2l Ul-y- )Kl-K4 

— S4 + m\ — m\, 



(36) 



A. Parametrizing the Cut Loop Momentum 



It is convenient to parametrize I in terms of the spinor rep- 
resentations of two massless real momenta which can be con- 
structed from any two adjacent external momenta on the box. 
(This parametrization was previously used in Refs. [60, 61, 
I65I1 .) I first choose two adjacent external momenta, K\ and 
K 2 and project them onto one another to form two massless 
momenta, K\ and k\- 



Ki= K\ 



Yn 



K 2 = K\ + 



Yn 



K 



1 ! 



K = 



J _ S i S 2 



(33) 



where 



Si — K\K\ 

Yn=2KlK b 2 



S 2 — K 2 -K 2l 
-K l -K 2 ±y/A(K h K 2 ), 



(34) 



A(K U K 2 ) = (Ki-K 2 f-SiS 2 . 



Note that there are two solutions for 712 unless either S\ = 
or S 2 = 0. 

Assuming that the Gram determinant, A(Ki,K 2 ), does not 
vanish, I can now form four massless vectors to use as a basis 
for solving for the loop momentum I. 



a 2 a^ 



3 2 — 
h OC4 ^4 ■ 



(35) 



ai a 2 



Yn 



a 3 = -- 



fo±\lPl-2P A TvlfL\]t^ 2 i 



2{K\-\Z A \K\- 



CLa, 



4a 3 



There are two solutions for 0:3 and it might appear that, com- 
bined with the two solutions for 712, there are four solutions 
for £. However, it works out that 



nriW)= ^(7r 2 ,« 3 -), 

Yu, 



(37) 



so that there are only two solutions for . 

This parametrization of V 1 looks rather different than that in 
Ref. 14111 . which used as a basis the three independent exter- 
nal momenta, , and and the antisymmetric com- 
bination of those three, P* 1 = e^ vpx K\ v K 2p K A x. Numer- 
ically, of course, the solutions are identical, but that given 
here has a number of features to recommend it. First, it is 
quite compact and easy to compute, even allowing for ar- 
bitrary values for the internal and external masses. Sec- 
ond, it is better behaved in Gram-singular configurations. As 
with all integral reduction techniques, we find spurious sin- 
gularities in the form of inverse powers of the Gram deter- 
minant arising from tensor reduction. The Gram determi- 
nants, A{K\ : K 2l Kn), are found in the on-shell solution to the 
loop mo mentum £. N umerical analysis shows that I scales 
like 1 / y/ 'A(K\,K 2 ,K4) near the Gram singularity. This prop- 
erty is explicit in Eq. i36\ , where the denominator of 0:3, 



K[ |^4 1 K\ j is the (complex) square root of A(K\ , K 2 , K4) . 

In the solution of Ref. Ir4lll . the components of I have coeffi- 
cients that scale like 1 /P 1 , where P 2 = A(K ly K 2 ,K4). Since I 
actually scales like 1 / A{K\ : K 2: Ka), there are cancellations 
hidden in the parametrization. This makes the problem of lo- 
cating and canceling these spurious singularities much more 
difficult. 



V. THREE PARTICLE CUTS AND THE EXTRACTION OF 
SCALAR TRIANGLE COEFFICIENTS 



This basis has the property that 03 -04 — 4a \ -a 2 — 2y\ 2 , but all 
other inner products of the a, among themselves vanish. The 



The procedure for extracting the box coefficients is so re- 
markably simple that one is inspired to try the approach for 



extracting the coefficients of lower point integrals. The use of 
triple cuts in one loop amplitudes to extract the triangle coef- 
ficient was first discussed by Mastrolia [66]. A complication 
arises when one imposes a triple cut on a one-loop amplitude. 
One does pick out a single triangle topology, but one also gets 
contributions from box topologies formed by splitting one of 
the triangle's vertices to open up a fourth propagator. Since 
the box contributions can be determined by the prescription 
above, one needs a way of separating the pure triangle con- 
tributions from the already known box terms. Forde B55I1 has 
recently described an elegant way of doing so. The discussion 
below follows that of Forde, but allows for massive internal 
legs on the triangles. 



The integral depicted in figure [3]is given by 

CjC (m ,-K 1 ,mi,K 2 ,m2) 

d 4 £ A[ n (Ki ; £)A ( 2 j) {K 2 ;£)Af {K 3 ; I) 
{2nf (P-mDde-K!) 2 - m\)((l + K 2 ) 2 - mj) 



n{-2ni) 



J ^A^\Kv,l)A'i>{K 1 -4)A^{Kr,i) 



iCJ) 



8(P-n%)8{{l-KiY-m\)8[{l + K 2 ) 2 -n%) 



i{-2%i) 



■J t A^{t)A 2 J> (t)A^{t) : 



(40) 




FIG. 3: a triple cut. 



I consider a triple-cut triangle, as shown in figure [3] The 
three delta function constraints imposed by the cuts, 

8(f-m 2 ), SaC-Ktf-m 2 ), 8((£ + K 2 ) 2 -m 2 2 ), (38) 



are not sufficient to completely fix the loop momentum, £, 
which must therefore have an unconstrained degree of free- 
dom. Using the same basis involving K\ and K 2 as for the 
box, I can parametrize the loop momentum as 



= ai af + a 2 a'i 



0=4 u 



OEi = 



0=2 



S2 (7i2 - Si ) + (yi2 - S 2 ) m 2 - 712 m\ + S 2 m\ 

H 2 -s l s 2 

S\{7n-S 2 ) + {yi 2 -Si)ml-y n m\ + Sim 2 2 
7 2 n ~SiS2 



(39) 



where J t is the Jacobian of the transformation from the delta- 
function constrained integral over d A i to the integral over the 
remaining free parameter t . Treating t as a complex variable 
and partial fractioning off terms with poles at finite t , this last 
integral can be rewritten as 



lnf t A[ j) A { 2 j) A { 3 i] ] (f) 
t-h 



(41) 



which represents a sum over the residues of all poles {k] at 
finite tic and a contribution at infinity. The [Inf, ] term is defined 
so that 

lim([Inf,A 1 A 2 A3](f)-Ai(f)A 2 (f)A 3 (f))-0. (42) 
The [Inf r ] term will be some polynomial in t , 



[Inf f AxA 2 A 3 ] (?) = 

h=0 



(43) 



where m is set by the maximum tensor rank allowed. For 
renormalizable theories, the maximum tensor rank for the tri- 
angle is three. 

Now, recall from the discussion of the box coefficient that 
the fourth delta function constraint fixed the value of t at some 
complex value to. Thus, the sum over residues at finite t sim- 
ply correspond to the contributions to the triple-cut from the 
various box configurations that satisfy those three cuts. More- 
over, the terms only give contributions to the scalar box co- 
efficients. The triangle contributions of the triple-cut come 
exclusively from the terms at infinity. That is, 



dt 



WtA^A^A^ 



(t) . (44) 



So now, I need to construct a table for integrals of the form 

dt 



-2%if 



-j t t n . 



(45) 



0=4 = a,\ a 2 



This is easily done by considering tensor triangle integrals of 
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the form 



d 4 £ 



(e-a 4 ) n 



{2K f (£2 _ m 2) ( (£ _ ^ )2 _ „2J ( (£ + tf 2 ) 2 - m 2 ) 

(-2*i) 3 ft ' 



(2*)' 



For the scalar integral, n = 0, it is clear that 



i{-2ni) 



dt 



J, = l. 



(46) 



(47) 



For tensor triangles however, with n > 0, the integrals vanish. 
By Lorentz invariance, the components of the tensor triangles 
must involve only products of K\, K2 or the metric tenor g^ v . 
But K\ and K2 decompose into K\ and K\, which annihilate 
via the Dirac equation, and metric tensor terms vanish because 
04-04 = 0. Therefore, 



and 



cj= - lnf t A[ j) A { 2 j) A { 3 j) 



C0U=o- 



(48) 



(49) 



VI. 



TWO PARTICLE CUTS AND THE EXTRACTION OF 
SCALAR BUBBLE COEFFICIENTS 



To extract the coefficients of bubble integrals, I proceed as 
before and impose the cuts that define the bubble topology 



8{f 



(50) 



Only one bubble configuration will satisfy these cuts, but mul- 
tiple triangle and box configurations will do so. Since the 
boxes and triangles can be extracted by the methods above, 
the task here is to isolate the pure bubble contributions. 




FIG. 4: a double cut. 

Since I only have one external momentum, K\, in a bubble 
configuration, I will choose an arbitrary massless momentum 
to define my parametrization. My massless projection of 
K\ and the set of basis vectors for the loop momentum is then 



defined in terms of 



Ki = K\ 4 



Si 



-x, 



y lx =2K l - X =2Klx, 



4 = z M - 



Kl\f\X- 



< = (x-\y"\k\ 



(51) 



: yai- 



- OC2 «2 



where both y and t are free parameters. The other coefficients, 
OC2 and 0C4 are fixed by the cut conditions and are found to be 



« 2 



S\+mk — m? S 



-y, a 4 



1 

2t 



Yix Tlx 

The integral depicted in figure |4]is 

biB (m ,-Ki,mi) 

d 4 £ Al\K 1 ;e)A^(K 2 ;e) 



a 2 y- 



(52) 



{2nf 
i{-2nif 



m 



t-^AfiK^AfiK^Dx 



{2% 



(53) 



= i (-23T if J J,, y Af (t,y)A® (t,y), 

where J ty is the Jacobian of the transformation. This time 
treating both y and t as complex variables and applying partial 
fractioning, this last integral can be rewritten as: 

bi = i{ -2 % if f^j t jAf{t,y)Af{t,y) 



(t,y) 
(y) 



(2nY 
dtdy 



InfvE 



E 

{./} 



{;} 



Res 



Resj^A^' A2 



t-t k 
\rf t AfAf 



y=yj 



(0 



y-yj 



(54) 



Res v=yj Z{k} 



(0„(0 



Res t= i k A\' A 



t-tk 



y-yj 



The double [Inf ] terms are the pure bubble contributions to the 
double-cut, the single residue terms are triangle contributions 
and the double residue terms are box contributions. As before, 
the box contributions only give the scalar box coefficients and 
are therefore not of interest to extracting bubble coefficients. 
The triangle terms are not so easily dismissed. Unlike in the 
previous section, the parametrization of the loop momentum 
does not annihilate all tensor triangle contributions. Since ten- 
sor triangles can be decomposed into scalar triangles, bubbles 
and tadpoles, there is a contribution to the bubble coefficient 
from the single residue terms. 
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A. Pure Bubble Contributions to the Bubble Coefficient 



Veltman reduction by using another set of identities, 



To extract the pure bubble contribution to the coefficient, I 
must build an integral table for powers of both t and y. If I 
consider tensor bubbles of the form 



/ 



{£-a 4 ) n 



{2nf {P-mD^l-Kyf-mX) 
dtdy 



(27Z) 



(55) 



d 4 £ 



{£ ■ a 3 ) k {£ ■ a 4 ) k 



(2tt) 4 {I 2 -ml) {{i-Kxf-mfj 



= {-\) k 2 k k\^B^ 



Si+ntQ — m\ 2 



2 k k\B^. 



(59) 



For n = 0, this is just the scalar bubble and the result is of 
course equal to unity. For n > 0, the integrals vanish because, 
by Lorentz invariance, tensor bubbles have components made 
up of Klf' and the metric tensor g^J, both of which annihilate 
products of ... . Thus, 



/ 



dtdy 
dtdy 

(27T 



Jt,y= 1: 



(56) 



The integral table for values of y is somewhat more compli- 
cated. First, I form the auxiliary vector 



Eliminating the Passarino-Veltman coefficient, I find 



dtdy 
dtdy 
dtdy 



Jt.y= 1. 



S\y 



2k- 1 



0, 



{2k) 



4 J ',y 



2 k k\ {Si+ml-m 2 
(2*j! 

1 



2Si 



2k 



2 k k\ 



m Q 



Si + OTq — m\ 
~Si 



= 0. 
(60) 



Tlx 

1 



Ki- (0 = 0, 



(Q- 



l-(0=^(Si+ml-mf)-Siy 



Using (0, 1 find, 



f d 4 £ 
J (2% 



(I- a) 



2k- 1 







-Si, 



(2k) 4 (£ 2 -m 2 ) ((l-Ktf-mft 

=> jdt dyJ t . y f j - Si y J = 

d 4 £ (£-(o) 2k _{2k)\ 



(57) 



(2k) 4 (£ 2 -m 2 Q ) ((l-Ktf-mft 2 k k\ 



(Sy) k B^ 



dtdyJ t .y 



S\ +ml — m 2 
2S^ 



Si J 2 k k\ 



(58) 



From these expressions, I can build an integral table for y: 



B(n) = 



dtdy 



{2k) 4 

(-1)" L '^ 2J 
n + l & Q 



Jt,y y" 



n — i 



—ml \' f S\ + ml — m\ 



n— 2i 



(61) 



B. Triangle Contributions to the Bubble Coefficient 



It is also possible for there to be triangle contributions to the 
bubble coefficient. These will come from the single residue 
terms in Eq. (l54b . The residue terms correspond to an addi- 
tional propagator going on-shell. Therefore, these terms can 
be obtained by applying an additional constraint to the double 
cut of the form, 



where B^° is the Passarino-Veltman reduction term contain- 
ing only products of the metric tensor (no factors of K\ ), and 
fi{° 2 *} i s the corresponding coefficient of the scalar bubble in 
that term. One can avoid actually constructing the Passarino- 



S{(£ + K 2 ) 2 -m 2 ) . 



(62) 



Without any loss of generality, y can be eliminated to satisfy 
the constraint, leaving t as the sole free parameter. The equa- 
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tion in y is 



C. The Bubble Coefficient 



= V 2 -^(,< K::\k; 



-ly 



K{K 2 --^ X .K 2 
l + niQ — m j 



2ty lx 



;rl* 2 |jq 



(63) 



The complete bubble coefficient is found by combining the 
pure bubble and triangle contributions. 



b, = -i 



Inf. 



Inf^A® 



(f,y) 



^0,y"^S(n) 



' triangles y—y± 



E E Inf f A«A«A« (0 



(67) 



t n ^T(n) 



where the A„ are the amplitudes formed by cutting one more 
propagator in either A y or A®. 



There are two solutions, y = y±, which must be averaged over. 
Because the parametrization of the loop parameter is based 
upon the vectors K\ and x, instead of K\ and K 2 , the tensor 
components of the triangle, that is, non-zero powers of t , will 
not identically vanish. 

The integral table for t can be worked out by considering 
tensor integrals of the form 



dH 



(l-a 4 ) n 



{!%)* - «§)((* -tfi) 2 - 
dtdy_ , 
(2nr J '< ■ 



>{-2Ki?f lx 



(64) 



Because 04 is a null vector and is orthogonal to K\, if I 
perform Passarino-Veltman reduction, I find that only the 
(K 2 ) component contributes to this integral. I then pick 
out the coefficient of the particular bubble I am looking at, 
Bo(mo,— K\,m\), from that Passarino-Veltman term. I find 
that 



T{n) = I (2^ % 



Sl y(x-\f2\Kl-) (K\ -K 2 ) 



(65) 



2y lx J A" (Ki,K 2 ) 

where A(Ki ,K 2 ) is the triangle Gram determinant and 



^1 

^2 



0, 
1, 
3 
2 



5 fS 



53 = 2 



A 

Si 

™2 



K\ K 2 



s 2 - 



Si 

2A{K U K 2 ) 
3 {K x -K 2 ) 2 



K\ K 2 
S\ + m 2 ) — m 



(66) 



VII. SINGLE PARTICLE CUTS AND THE EXTRACTION 
OF SCALAR TADPOLE COEFFICIENTS 

The tadpole coefficients can be extracted by an extension of 
the same procedure. The only constraint that the loop momen- 
tum must satisfy is that 1 2 — m 2 , = 0. Satisfying this constraint 
leaves three free parameters in £, 



at, 



It 



1 / m. 
wy - 



(68) 



where %^ and xj/f 1 are arbitrary light-like momenta, since, in 
a tadpole configuration, there are no external momenta on 
which to base the parametrization. 




i- 1 



FIG. 5: a single cut. 
The integral depicted in figure [5] is 



a m Ao(mo) = i 



(2n) 4 (£ 2 -mh 



i{-2%i) 



{in) 



r A( m )(0;£)5(^ 2 -m 2 ) (69) 



where A(0;€) indicates that there is no external momentum 
flowing in to the tadpole. I again treat the free parameters as 
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. f dtdydw 
( >J ~EnW Jt ' y > w: 



complex variables and apply partial fractioning to obtain, 

( () J (2nf • / "-- A ' : '-- v -"' : ' = 

"dtdydw (70) 

([Inf w [Infy [Inf, A]]] w) + ...). 

The triple [Inf ] term will be the pure tadpole contribution to 
the tadpole, the double [Inf] - single [Res] terms will be the 
bubble contribution to the tadpole and the single [Inf ] - double 
[Res] terms will be the triangle contribution to the tadpole. 
As before, the pure residue terms correspond to scalar box 
contributions and need not be considered. 



A. Pure Tadpole Contributions to the Tadpole 

To extract the pure tadpole contribution to the coefficient, 
I must build an integral table for powers of w, y, and t. This 
task is greatly simplified by the fact that tensor tadpole inte- 
grals vanish when the rank of the tensor is odd. Since there are 
no external momenta to project upon, the only non-vanishing 
tensors are proportional to products of the metric tensor. Fur- 
ther simplifications follow from the fact that the basis vectors 
are light-like and the only non-vanishing contractions are 

(x~\y p \w~)(w~\yh\x~) = -^x-w- ( 71 ) 

These facts lead to the conclusion that the only non-vanishing 
entries in the integral table are of the form 

' dtdydw 



(27T) 



■J t , y , w y n w n ^0, 



(72) 



where n > 0. All integrals involving non-zero powers of t 
vanish identically. 



/ 



/ 



d 4 £ (£-x) k (£-W) k 
(2nf £ 2 -ml 

-<-2,o (to) 



a(M) A w> 



2k 



d 4 t Or|/|vT(v 



dtdydw 

\x-) k 



Jt,y,w (wy) k 



(2Kf 



£ 2 - 



kl (-2 7zv ) k AW 



dtdydw 

l2^Y~ Jt ' y - w 



wy - 



2 \ * 

mi) 

(73) 



Together, these equations yield the result that 

dtdydw 



(_2,r '" ) / : W 1/ '' y,w= 1 

(-2*0 j dtdydW 



(2nY 



-J, 



(-wy) 1 



A consistent solution to this set of equations is 



D(k) = j 



dtdydw 
l2^F 



Jt, y ,w {wy) k = 



1 



k+l 



m Q 

Yxv 



0. 

(74) 

(75) 



B. Bubble Contributions to the Tadpole 

To obtain the bubble contributions to the tadpole, I add a 
second constraint, 

5((* + *i) 2 -m?). (76) 
This constraint can be imposed by solving 

0=y W i^ ^'1* I + y 2x-Ki+w2y-K l 

+ t( X -\^\ W -)-^( W -\^\ X -) W 



tYxv 



+ S\+m Q -m l 



for the values w = wq or y = yo and then constructing an in- 
tegral table for t and y or t and w, respectively. In fact, one 
must average over the two solutions. The integral tables are 
constructed by considering tensor integrals of the form 



/ 



d 4 £ 



(£-y) k (y- 



(2nf (P-ml)((£ + K x ) 2 -m 



W=WQ 



r dtdy_ k 

4 J t.yJ 1 



I 



d 4 £ 



>{-2*if(J-f) (- r „) 
(£-x) k (v-\i\x-) m 



(2*)< 



(78) 



(2tt) 4 {£ 2 -ml){(£ + Kyf-m^ 

= (K 1 . X ) k ( ¥ -\K 1 \ X -) m B {lk+a} 

-(-2n l f(^)\-rx V rj^A^ 

The integral tables are: 

Mk,n-k) = j^J' t , y y k ?- k = 

(-If 1 f2 W - Kl \ k f (yr\Zi\X-y 



y=yo 



n + l Si V 7x¥ 
L(--iV 2 J /„-l-i 
i 



i=0 



Si 



Txv 
Si +mQ — m\ 



n-\-2i 



y (k,n-k) = J 



(2nf 
(-1)" 1 (2%-Ki 



dtdy i k n -k _ 



« + l Si V Yxv 
L(»-i)/2j /-_!_.• 

s ( . 



(W-\Ki\x-) 



Yx¥ 



n—k 



Si 



Si + m^ — m\ 

s7 



2 \ n-1-2/ 



(79) 



C. Triangle Contributions to the Tadpole 

The triangle contributions to the tadpole come from adding 
a third constraint, 



8((£ + K 2 ) 2 -ml) 



(80) 
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The extra constraints can be imposed by simultaneously solv- 
ing Eq. ( |77] i and 

= ywi^-iM^J. +y 2% ■ K 2 + w2\j/ ■ K 2 



D. The Tadpole Coefficient 

The complete tadpole coefficient is found by combining the 
pure tadpole, bubble and triangle contributions. 



■t{x-\ti\vr)-^-{vr \*%\z~) 

r rx¥ 

- S 2 + niQ — m 2 



(81) 



for y and w. There are two solutions, (w,y) = (w\,yi) and 
(w,y) = (w 2 ,y 2 ), which must be averaged over. This leaves 
one free parameter t , for which I must derive an integral table. 
This can be done by considering tensor integrals of the form 

Or War)" 

{2nf (*-m§) ((l + Ktf-ml) ((£ + K 2 ) 2 - m 2 ) 



i=0 



(2; 

where C^ 1 " ' 2 '^ is the n-th rank Passarino-Veltman triangle 
coefficient whose Lorentz structure has (n — 1) powers of K\ 
and i powers of K 2 . I only need the tadpole term from the 
Passarino-Veltman coefficients, so my integral table is: 



F(n) 



dt 



■ f t 't n 



" {y-\Kr\x-)"-{y-\K2\X-) 1 f n 
i=0 1 xv 



(83) 



with the / n (()'s given by 

/o(0) =0, 
/i(0) =/i(l)=0 
1 K x -K 2 



fi(0) 
/ 2 (2) 
/ 3 (0) 



/2(1) = - 



1 



1 



2A(K U K 2 ) 



4 SiA(Ki,K 2 ) ' 
1 K x -K 2 
~~4S 2 A(K h K 2 ) ' 

1 {Sx+ml-m^jKyK^ 
6 S 2 A{K U K 2 ) 
\S 2 + m 2 ) — m 2 l 5 S 2 (S 2 + m 2 ] — mi ! ) 
8 S X A(K U K 2 ) + 24 A 2 {K U K 2 ) 
5 S 2 (S 1 +m 2 -m 2 1 )(K 1 -K 2 )^ i84) 



/ 3 (1) = 



SiA 2 (K u K 2 ) 

0~"'l 



24 

/ 1 Si + m 2 ) — m 2 
\24 S\A(Ki,K 2 ) 

5 {S 2 + m 2 -m 2 2 )(Ki-K 2 ) 
24 A 2 (K U K 2 ) 
5 S 2 (Si+m 2 ~m 2 ) 
24 A 2 (K U K 2 ) 

/3(2) =/3(l)| J f 1 ^ 2 , mi „ OT2 , 
/3(3) =/3(0)|^^ 2 !B|H „ 2 ■ 



Inf„, 



Inf v 



Inf,A ( 



(t,y,w) 



bubbles 



Inf v 



Infw 



Inf f A] 'A\ ' 



f-»0,(wy)"-»D(«) 



y k t'"^E w (k,m) 



lnf t A^lf 5 



y=vo 



(w 2 ,)'2) 



' triangles (ny;) = (wi,;yi) 



where the A„ are the amplitudes formed by cutting a prop- 
agator in AW and A„ are the amplitudes formed by cutting 
two propagators in A . 



t"^>F{n) 



(85) 



VIII. COMMENTS 

The expressions given above assume that all external mo- 
menta are massive and that all internal masses are distinct. 
They have been tested numerically by decomposing tensor in- 
tegrals (through the fourth rank tensor box) using Passarino- 
Veltman techniques and those described here. 

As long as the external momenta are massive, the only com- 
plication in allowing the internal masses to become degenerate 
comes from the tadpoles. In case of degeneracy, one must be 
careful not to over count the contributions from bubbles and 
triangles. This can clearly happen because the starting ampli- 
tude for the tadpole cut A^"\%,i) is not tied to any external 
momenta. Thus, if one successively adds cuts to the tadpole 
configuration, one can reach the same bubble or triangle con- 
figuration from different starting points. It is safer to construct 
special cases for degenerate mass bubbles and triangles and 
simply sum over the distinct configurations. 

For S\ 7^ and niQ = m\ 7^ 0, 



(-1)" 1 f 2 W -KA k f (vr|j:, x 



7xv 



n—k 



L«/2J 

x I 

i=] 



n — 1 \ 1 —m, 



2\ ' 




Ey(k,n-k)= I ^yiyr^ 



(86) 



(-1)" 1 (2x-K x 



{2nf ,M ' 

k 



n+l Si V 7xv 

L«/2J 

x L 



(rj£ilr 

Yxv 



n-k 



n — i\ 1 —m, 
i ) [ a. 



2\ ' 
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Some formulae must also be modified when external mo- 
menta are massless. In particular, the bubble integral is no 
longer an independent loop-integral function, 



A (m ) -A (mi) 



Bo(0,mo,^i,mi)| Sl=0imo=mi ^ = (l-e 
flo(0,»io,£i,mi)| Sl=0 , mo=m ,=o = 0- 



o _m i 
A (m ) 



(87) 



In addition to obviating the need to extract the coefficient of 
the massless bubbles, this identity alters the extraction of tad- 
pole terms from bubbles and triangles. The expressions in 
Eq. j79l ) become for Si = , mo ^ m\, 



E w (k,n~k) = 



(-1)* 



n+l 



Ey(k,n — k) = 



{-If 



(2 V -Ki) k ( V -\l[i\x 



— \n—k 



(88) 



(2 X -Ki) k ( ¥ - \Kx\x~ 



\ n—k 



m oYi 



'XV 



and for Si = , mo = mi ^= 0, 

(-1)* (2vr-*i)*- 



E w (k,n 
Ey(k,n 



-k) = 



¥' 



-iir 



M —\n—k 



i 



-*) = 



(-1)* p*.^^- 



\n—k 



(89) 



1 



'XV 



Clearly this last case is one in which one must be careful to 
count the contribution of this bubble only once, since it can be 
reached two different starting points. 

Formulae which depend upon Passarino-Veltman triangle 
coefficients also change significantly. Actually, Eq. (l66l > is 
valid in the S2 — * limit and does not appear in the Si — > 
limit since, the bubble function So(mo, — ^i,mi)| s =0 is not 
a basis integral. Eq. d84i l. however, changes dramatically as 
the external momenta vanish and the internal masses become 
degenerate. There are five different special cases to consider 
as one takes different combinations of S„ = and m„ = 
(taking into account exchange symmetries). The full set of 
expressions is given in Appendix (lAb. 



The method is valid for arbitrary internal and external masses 
and can be applied to any one-loop calculation in quantum 
field theory. When combined with one-loop on-shell recur- 
sion relations, this procedure can be used to construct com- 
plete one-loop amplitudes. 

One of the features of this method is that it can make use of 
the compact representation of tree-level amplitudes found by 
using helicity methods. This formalism, when combined with 
on-shell recursion relations to determine the rational terms, 
as envisioned in the unitarity bootstrap, shows great promise 
for constructing an automated system for generating one-loop 
amplitudes. 



IX. CONCLUSIONS 



I have described a method, based upon generalized unitar- 
ity, for computing the four-dimensional coefficients of scalar 
loop-integral functions in one-loop amplitudes in quantum 
field theory. The result is that the coefficients of the loop 
integrals are determined from the product of tree-level dia- 
grams that appear at the vertices of the loop-integral itself. 
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APPENDIX A: SPECIAL CASES FOR THE TRIANGLE 
CONTRIBUTION TO THE TADPOLE 

1. Case 1: Si = 0, m\ ^ mo 

Because the massless scalar bubble breaks into a sum of 
tadpoles, the special cases gives non-vanishing contributions 
toF(l). 



/o(0) 
/l(0) 

/ 2 (0) 



/ 2 (1) 

Mo) 



2{K l -K 2 )(m 2 -m 2 ) 



H) 



4 \{K l -K 2 ){ml-m\) 2 
3 S 2 



/l(l)=0 

s 2 + 



(K l -K 2 ) 2 (m 2 - 



2{K l -K 2 f 



%{Ki-K 2 ) 2 
1 1 
"6 {Ky-K^m 2 
1 (S 2 + m 2 - m 2 ) 2 25 S 2 {S 2 + m 2 - m 2 ) 



1 {S 2 + ml-m\) 
8 m 2 JK x -K 2 ) 2 



8 m 2 (Kf K 2 f 
5 S\(m\-m 2 ) 
16 (K\ -K 2 ) 5 ' 



48 {Ki-K^ 4 



Mi) = 

/3(2) = 
/3(3) 



5 S 2 (m 2 — m 2 ) 
16 (K x -K 2 f 
1 



1^2- 

3 {Kr-Ktf 
n 2 



S(m 2 -mi)(K l -K 2 y 



1 mi 



1 S 2 + ml- 



4(K!-K 2 ) 3 24 S 2 (K l -K 2 ) 2 



1 (m 



1S 2 - 



6 £{K X -K 2 ) 



(Al) 



2. Case 2: Sj = 0, wi = mo 

This is perhaps the most difficult case and can be tricky, 
since the invariant mass of the third external momentum, 
K 3 = K\ — K 2 appears in the denominator. It can happen that 
53 = 0, in which case one should be looking at Case 4 (see Ap- 
pendix (IA41 i). and if one had arrived at this configuration by 
cutting the other leg of the triangle first, one would naturally 
have arrived there. Therefore, the formula below assumes that 
S 3 = (K l -K 2 ) 2 ^0. 



/o(0) 
/i(0) 

/ 2 (o) 

/ 2 (1) 

/ 2 (2) 
/ 3 (0) 



=0, 
1 



2m\ {K x -K 2 ) 

=0, 



1 



1 



1 (K 1 -K 2 ) + 2m, 
4 



- S 2 - mi 



4S 3 (K l -K 2 ) 
1 1 
"4Si{Ki-K 2 y 
1 1 



AS 2 {K X -K 2 ) 4 S 3 (K\ -K 2 ) ' 
1 7 53 — 2 (otq — m 2 ) 1 Wq — 



2 S 3 {K X -K 2 f 



8S3 (K v K 2 y 



1 4 ( Ki ■ K 2 ) + 1 4 ml + 3 S 2 - 3 m\ 
24 m 2 (K x -K 2 ) 2 
1 {S 2 -m 2 ) 2 1 l9S 2 + Sml-Um 2 , 



8m2(^.^ 2 ) 3 24 {K x -K 2 f 



/3(1) 
/ 3 (2) 

/ 3 (3) 



1 2S 3 



1 5S 3 - 



6 Sf(^i.AT 2 ) 
1 S3 — 2 (otq — mj) 
12 S 2 {KyK 2 ) 



12 5 3 (^i-^2) 2 



24 5 3 (^i-^ 2 ) 2 24 5 2 (^i-^ 2 ) 2 
1 S3 + mo — ffJ 2 1 5 2 + iftk — m\ 



6 S 2 {K X -K 2 ) 6 S 2 {K X -K 2 ) 



(A2) 
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3. Case 3: S\ — S 2 — 0, mi 7^ rriQ , m 2 7^ mo 

If both Si = and S2 = 0, the third external momentum 
Kt, = K\ — K2 must be massive, S3 = — 2K\ ■ K2 7^ 0. In the 
limit that S3 becomes small, the diagram describes a one-loop 
splitting amplitude, rather than a scattering amplitude. 

/o(0) =0, 

MO) 1 
/i(i) 

Mo) 
Mi) 

M2) 

Mo) 



~2(m 2 -m 2 )(K 1 -K 2 ) ' 

1 1 
'2(m 2 -m 2 )(K l -K 2 ) ' 
1 ml 



1 



4 {m 2 -m 2 ) 2 (K l -K 2 ) 4 { m 2 - m 2 ){K, ■ K 2 ) 2 
1 1 



2(^1-^)2' 

1 ml 



1 



4(mg-ml) 2 (^i-«'2) 4(m 2 -m 2 )(tf 1 -tf 2 ) 2 

1 m o 1 m 2 {m 2 -m 2 ) 

"6 (ml-m 2 ?^ ■ K 2 ) 8 (m 2 -™ 2 ) 2 ^ ./T 2 ) 2 
1 (mg-m 2 / 



1 



/3(2) 



/3(3)=T 



(m 2 


-m 2 )(/Ti 


•^ 2 ) 3 




m 2 










(m 2 


-m 2 )(^i 






m 2 





3 m 2 - 



3 niQ — mj 



S(m 2 -m 2 )(K l -K 2 ) 2 8 • K 2 ) 3 ' 

wiq 1 OTq(/Wq — m 2 ) 



e^-m 2 ) 3 ^-^) S^-m 2 ) 2 ^-^) 2 
(m 2 -m 2 ) 2 



1 

+ 77 



8(m 2 -m 2 )(^-/:2) 3 



(A3) 



4. Case 4: Si = S 2 = 0, mi = mo , m 2 7^ mo 
Again, in a scattering amplitude, S3 = — 2/q • /f 2 ^ 



/o(0) 
/i(0) 

/i(l) 

MO) 
Mi) 
M2) 
MO) 

Mi) 
M2) 

/3(3) 



1 



=0, 
1 

~2m 2 (K 1 -K 2 ) 
1 1 



^(mg-m 2 )^-^)' 
1 (K 1 -K 2 )-m 2 3 



1 



4 m 2 (K x -K 2 ) 2 i(K x -K2) 2 ' 

1 1 

4>i^ 2 ) 2 ' 



4(m 2 -m 2 ) 2 (K 1 -K 2 ) 8 (*i • tf 2 ) 2 ' 



1 



1 



■ + ; 



1 



6 m 2 (tfi.tfz) 24(^ r 

1 OTj (C^l ' ^2) + «o — m 2 ) _|_ 1 1 Wq — OT 



m§(*i-* 2 ) 3 



48 (Zi -^ 2 ) 3 



12(^1-^2) 



8 (^i-^ 2 ) 3 



1 



1 m 2 - 



S(m 2 -m 2 2 )(K 1 -K 2 ) 2 16 (K r K 2 ) 3 
m\ 1 2(^i ■K 2 )-ml + ml 



~6(m 2 -m 2 y 24 (^i-tf 2 ) 3 



(A4) 



5. Case 5: S\ = S2 — 0, m 2 = mi = mo 



If both Si = and S2 = 0, and all the internal masses are 
degenerate, 



/o(0) -0, 
/i(0) =/i(l) 

M0) =/ 2 (2) 

/ 2 (i)=o, 

/ 3 (0) =/ 3 (3) 

Mi)=M2) 



1 



1 



2m 2 (/: 1 -/: 2 ) 
1 1 



Am 2 (K l -K 1 ) 



+ ■ 



(A5) 



6 m 2 (K l -K 2 ) l2(Ki-K 2 ) 



2 - 



1 



1 



2 ■ 
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6. Case6: S\ =0, S2,S^ ^=0, ni2 = mi = niQ 7. Case 7: Si, 52, 53 ^= 0, m 2 = m\ = mo 

If only Si=0, but all the internal masses are degenerate, If ^ extemal masses are non . zero but all ±e internal 

, . masses are degenerate, 
7o(UJ =0, 

/'(°) = L.2/J FT ' /o(0)=0, 



/i(l)=0, 



/i(0)=/i(l)=0, 
/2(0)=/ 2 (l)=/ 2 (2)=0, 



1 1 1 S 2 

/2(0) - - - 2 -7, /jr . Kn \l ' A/^_ 15 2-(^f^2) 



4j^(^ 1 .AT 2 ) 4m 2 (K l -K 2 ^' / 3 ( ) 
/ 2 (l)=/ 2 (2)=0, 



6 S3A(^i,Ar 2 ) 

1 s 2 



^=64(^) + 6^(^) H^iT 1 '^ 65l( " 1 " 2) ' (A7) 

6 8 »8 i5 2 -i 3 -(*-*2) 

I 1 | 1 ^ /3(2) -6 S 3 A( M ) ' 

3(*i-Ar 2 ) 3 Smg^.^' l S 2 S3-^ + (g2 + S3)(gi-g 2 ) 

f,m=-I l - - l - 6 s 2 S3A(ati,a:2) 

731 ' 6S 3 {K X -K 2 ) 3(K l -K 2 )2> 
1 1 
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